Magnetic field decay of three interlocked flux rings with zero linking number 
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The resistive decay of chains of three interlocked magnetic flux rings is considered. Depending on the relative 
orientation of the magnetic field in the three rings, the late-time decay can be either fast or slow. Thus, the 
qualitative degree of tangledness is less important than the actual value of the linking number or, equivalently, 
the net magnetic helicity. Our results do not suggest that invariants of higher order than that of the magnetic 
helicity need to be considered to characterize the decay of the field. 
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I. INTRODUCTION 

Magnetic helicity plays an important role in plasma physics 
iil-lsfLsolar physics 0-0], cosmology 0-0], and dynamo the- 
ory Hit]. This is connected with the fact that magnetic 
helicity is a conserved quantity in ideal magnetohydrodynam- 
ics |[l2n . The conservation law of magnetic helicity is ulti- 
mately responsible for inverse cascade behavior that can be 
relevant for spreading primordial magnetic field over large 
length scales. It is also likely the reason why the magnetic 
fields of many astrophysical bodies have length scales that are 
larger than those of the turbulent motions responsible for driv- 
ing these fields. In the presence of finite magnetic diffusivity, 
the magnetic helicity can only change on a resistive time scale. 
Of course, astrophysical bodies are open, so magnetic helicity 
can change by magnetic helicity fluxes out of or into the do- 
main of interest. However, such cases will not be considered 
in the present paper. 

In a closed or periodic domain without external energy sup- 
ply, the decay of a magnetic field depends critically on the 
value of the magnetic helicity. This is best seen by consid- 
ering spectra of magnetic energy and magnetic helicity. The 
magnetic energy spectrum M{k) is normalized such that 



/ 



M{k)dk = (S^)/2/io, 



(1) 



where B is the magnetic field, /io is the magnetic permeabil- 
ity, and k is the wave number (ranging from to oo). The 
magnetic helicity spectrum H{k) is normalized such that 



/ 



H{k)dk = {A-B), 



(2) 



where A is the magnetic vector potential with S = V x A. 
In a closed or periodic domain, H{k) is gauge-invariant, i.e. 
it does not change after adding a gradient term to A. For fi- 
nite magnetic helicity, the magnetic energy spectrum is bound 
from below |[T^ such that 

M{k) > k\H{k)\/2fio- (3) 

This relation is also known as the realizability condition 
Thus, the decay of a magnetic field is subject to a correspond- 
ing decay of its associated magnetic helicity. Given that in a 



closed or periodic domain the magnetic helicity changes only 
on resistive time scales |14], the decay of magnetic energy 
is slowed down correspondingly. More detailed statements 
can be made about the decay of turbulent magnetic fields, 
where the energy decays in a power-law fashion proportional 
to t~^. In the absence of magnetic helicity, {A - B) = 0, we 
have a relatively rapid decay with a ^ 1.3 [15], while with 
(A • ^ 0, the decay is slower with a between 1/2 [9] and 
2/3 iPj. 

The fact that the decay is slowed down in the helical case is 
easily explained in terms of the topological interpretation of 
magnetic helicity. It is well known that the magnetic helicity 
can be expressed in terms of the linking number n of discrete 
magnetic flux ropes via 1.13.1 
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A-BdV = 2n$i$2, 



where 



is. 



B ■ dS (for z = 1 and 2) 



(4) 



(5) 



are the magnetic fluxes of the two ropes with cross-sectional 
areas Si and 5*2 . The slowing down of the decay is then plau- 
sibly explained by the fact that a decay of magnetic energy is 
connected with a decay of magnetic helicity via the realizabil- 
ity condition ([3]). Thus, a decay of magnetic helicity can be 
achieved either by a decay of the magnetic flux or by mag- 
netic reconnection. Magnetic flux can decay through anni- 
hilation with oppositely oriented flux. Reconnection on the 
other hand reflects a change in the topological connectivity, as 
demonstrated in detail in Ref. 1 17, p. 28]. 

The situation becomes more interesting when we consider 
a flux configuration that is interlocked, but with zero linking 
number. This can be realized quite easily by considering a 
configuration of two interlocked flux rings where a third flux 
ring is connected with one of the other two rings such that the 
total linking number becomes either or 2, depending on the 
relative orientation of the additional ring, as is illustrated in 
Fig.[T] Topologically, the configuration with linking numbers 
of and 2 are the same except that the orientation of the field 
lines in the upper ring is reversed. Nevertheless, the simple 
topological interpretation becomes problematic in the case of 
zero linking number, because then also the magnetic helicity 



2 



form 






(6) 
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FIG. 1: Visualization of the triple ring configuration at the initial 
time. Arrows indicate the direction of the field lines in the rings, cor- 
responding to a configuration with n = (left) and n = 2 (center). 
On the right the non-interlocked configuration with n = is shown. 



is zero, so the bound of M from below disappears, and M can 
now in principle freely decay to zero. One might expect that 
the topology should then still be preserved, and that the link- 
ing number as defined above, which is a quadratic invariant, 
should be replaced with a higher order invariant 1 18-20]. It 
is also possible that in a topologically interlocked configura- 
tion with zero linking number the magnetic helicity spectrum 
H{k) is still finite and that bound ^ may still be meaning- 
ful. In order to address these questions we perform numerical 
simulations of the resistive magnetohydrodynamic equations 
using simple interlocked flux configurations as initial condi- 
tions. We also perform a control run with a non-interlocked 
configuration and zero helicity in order to compare the mag- 
netic energy decay with the interlocked case. 

Magnetic helicity evolution is independent of the equation 
of state and applies hence to both compressible and incom- 
pressible cases. In agreement with earlier work [21] we as- 
sume an isothermal gas, where pressure is proportional to den- 
sity and the sound speed is constant. However, in all cases the 
bulk motions stay subsonic, so for all practical purposes our 
calculations can be considered nearly incompressible, which 
would be an alternative assumption that is commonly made 



~Dt 



Dlnp 



(8) 



where Fvisc = p • 2vpS is the viscous force, S is 
the traceless rate of strain tensor, with components Sij = 



V X B/po is the current 



density, v is the kinematic viscosity, and r] is the magnetic 
diffusivity. 

The initial magnetic field is given by a suitable ar- 
rangement of magnetic flux ropes, as already illustrated 
in Fig. [TJ These ropes have a smooth Gaussian cross- 
sectional profile that can easily be implemented in terms 
of the magnetic vector potential. We use the Pencil Code 
(http : / /pencil-code . googlecode . com), where 
this initial condition for A is already prepared, except that 
now we adopt a configuration consisting of three interlocked 
flux rings (Fig. [T]) where the linking number can be chosen 
to be either or 2, depending only on the field orientation in 
the last (or the first) of the three rings. Here, the two outer 
rings have radii Rq, while the inner ring is slightly bigger and 
has the radius Ri = 1.2Ro, but with the same flux. We use 
Ro as our unit of length. The sound travel time is given by 
Ts = Ro/cs. 

In the initial state we have U = and p = po = I. Our 
initial flux, ^ = J B dS, is the same for all tubes with 
^ = OA CsR^^/poPO' This is small enough for compressibil- 
ity effects to be unimportant, so the subsequent time evolu- 
tion is not strongly affected by this choice. For this reason, 
the Alfven time, Ta = ^^poPoRl/^, will be used as our time 
unit. In all our cases we have Ta = lOTg and denote the di- 
mensionless time as r = t/TA- In all cases we assume that 
the magnetic Prandtl number, u/rj, is unity, and we choose 
u = T] = IQ-'^RoCs = 10-^RI/Ta. We use 256^ mesh 
points. 

We have chosen a fully compressible code, because it is 
readily available to us. Alternatively, as discussed at the end of 
§111 one could have chosen an incompressible code by ignoring 
the continuity equation and computing the pressure such that 
V • t/ = at all times. Such an operation breaks the locality 
of the physics and is computationally more intensive, because 
it requires global communication. 



II. MODEL 



III. RESULTS 



We perform simulations of the resistive magnetohydrody- 
namic equations for a compressible isothermal gas where the 
pressure is given by p = pc^, with p being the density and 
Cg being the isothermal sound speed. We solve the equations 
for A, the velocity Z7, and the logarithmic density In p in the 



Let us first discuss the visual appearance of the three inter- 
locked flux rings at different times. In Fig. [2] we compare the 
three rings for the zero and finite magnetic helicity cases at 
the initial time and at r = 0.5. Note that each ring shrinks as 
a result of the tension force. This effect is strongest in the core 
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FIG. 2: Visualization of the triple ring configuration at r = (left), 
as well as at T = 0.5 with zero linking number (center) and finite 
linking number (right). The three images are in the same scale. The 
change in the direction of the field in the upper ring gives rise to a 
corresponding change in the value of the magnetic helicity. In the 
center we can see the emergence of a new flux ring encompassing 
the two outer rings. Such a ring is not seen on the right. 



of each ring, causing the rings to show a characteristic inden- 
tation that was also seen in earlier inviscid and non-resistive 
simulations of two interlocked flux rings |21]. 

At early times, visualizations of the field show little differ- 
ence, but at time r = 0.5 some differences emerge in that the 
configuration with zero linking number develops an outer ring 
encompassing the two rings that are connected via the inner 
ring; see Fig.[2l This outer ring is absent in the configuration 
with finite linking number. 

The change in topology becomes somewhat clearer if we 
plot the magnetic field lines (see Fig. O. For the n = 2 con- 
figuration, at time r = 4 one can still see a structure of three 
interlocked rings, while for the n = case no clear structure 
can be recognized. Note that the magnitude of the magnetic 
field has diminished more strongly forn = than forn = 2. 
This is in accordance with our initial expectations. 

The differences between the two configurations become 
harder to interpret at later times. Therefore we compare in 
Fig. [4] cross-sections of the magnetic field for the two cases. 
The xy cross-sections show clearly the development of the 
new outer ring in the zero linking number configuration. From 
this figure it is also evident that the zero linking number case 
suffers more rapid decay because of the now anti-aligned 
magnetic fields (in the upper panel is of opposite sign 
about the plane y = while it is negative in the lower panel). 

The evolution of magnetic energy is shown in Fig.[5]for the 
cases with zero and finite linking numbers. Even at the time 
r ^ 0.6, when the rings have just come into mutual contact, 
there is no clear difference in the decay for the two cases. 
Indeed, until the time r 2 the magnetic energy evolves 
still similarly in the two cases, but then there is a pronounced 
difference where the energy in the zero linking number case 
shows a rapid decline (approximately like t~^/^), while in the 
case with finite linking number it declines much more slowly 
(approximately like t~^/^). However, power law behavior is 
only expected under turbulent conditions and not for the rel- 
atively structured field configurations considered here. The 




FIG. 3: (Color online) Magnetic flux tubes at time r = 4 for the 
case of zero linking number (upper picture) and finite linking num- 
ber (lower picture). The colors represent the magnitude of the mag- 
netic field, where the scale goes from red (lowest) over green to blue 
(highest). 



energy decay in the zero linking number case is roughly the 
same as in a case of three flux rings that are not interlocked. 
The result of a corresponding control run is shown as a dot- 
ted line in Fig. [S] At intermediate times, 0.5 < r < 5, the 
magnetic energy of the control run has diminished somewhat 
faster than in the interlocked case with n = 0. It is possible 
that this is connected with the interlocked nature of the flux 
rings in one of the cases. Alternatively, this might reflect the 
presence of rather different dynamics in the non-interlocked 
case, which seems to be strongly controlled by oscillations on 
the Alfven time scale. Nevertheless, at later times the decay 
laws are roughly the same for non-interlocked and interlocked 
non-helical cases. 
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FIG. 4: (Color online) Cross-sections in the xy plane of the magnetic field with zero linking number (upper row) and finite linking number 
(lower row). The z component (pointing out of the plane) is shown together with vectors of the field in the plane. Light (yellow) shades 
indicate positive values and dark (blue) shades indicate negative values. Intermediate (red) shades indicate zero value. 



The time when the rings come into mutual contact is 
marked by a maximum in the kinetic energy at r 0.6. This 
can be seen from Fig.[6l where we compare kinetic and mag- 
netic energies separately for the cases with finite and zero link- 
ing numbers. Note also that in the zero-linking number case 
magnetic and kinetic energies are nearly equal and decay in 
the same fashion. 

Next we consider the evolution of magnetic helicity in 
Fig. [71 Until the time r ^ 0.6 the value of the magnetic helic- 
ity has hardly changed at all. After that time there is a gradual 
decline, but it is slower than the decline of magnetic energy. 
Indeed, the ratio {A-B) / {B'^), which corresponds to a length 
scale, shows a gradual increase from 0.1i?o to nearly O.QRq at 
the end of the simulation. This reflects the fact that the field 



has become smoother and more space-filling with time. 

Given that the magnetic helicity decays only rather slowly, 
one must expect that the fluxes of the three rings also only 
change very little. Except for simple configurations where 
flux tubes are embedded in field-free regions, it is in general 
difficult to measure the actual fluxes, as defined in Eq. ([5]). 
On the other hand, especially in observational solar physics, 
one often uses the so-called unsigned flux lHHIl^], which is 
defined as 



2D 



Bids. 



(9) 



For a ring of flux ^ that intersects the surface in the middle at 
right angles the net flux cancels to zero, but the unsigned flux 
gets contributions from both intersection, so P2D = 2 1 $ | . In 
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FIG. 5: Decay of magnetic energy (normalized to the initial value) 
for linking numbers 2 (solid line) and (dashed line). The dotted 
line gives the decay for a control run with non-interlocked rings. The 
dash-dotted lines indicate t^^^ and t^^^ scalings for comparison. The 
inset shows the evolution of the maximum field strength in units of 
the thermal equipartition value, Bth — Cs(poMo)^^^- 
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FIG. 6: Comparison of the evolution of kinetic and magnetic en- 
ergies in the cases with finite and with vanishing linking numbers. 
Note that in both cases the maximum kinetic energy is reached at the 
time T ^ 0.6. The two cases begin to depart from each other after 
T 2. In the non-helical case the magnetic energy shows a sharp 
drop and reaches equipartition with the kinetic energy, while in the 
helical case the magnetic energy stays always above the equipartition 
value. 



three-dimensional simulations it is convenient to determine 

P= j \B\dV. (10) 
Jv 



FIG. 7: Evolution of magnetic helicity in the case with finite linking 
number. In the upper panel, {A - B) is normalized to its initial value 
(indicated by subscript 0) while in the lower panel it is normalized to 
the magnetic energy divided by Ro. 
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FIG. 8: Decay of the unsigned magnetic flux P (normalized to the 
initial value Pq) for the cases with n = and n = 2. The dotted line 
gives the decay for a control run with non-interlocked rings. 



For several rings, all with radius R, we have 

N 

P = 27rR^\^i\=7rNRP2B. (11) 

i=l 

where N is the number of rings. In Fig.[8]we compare the evo- 
lution of P (normalized to the initial value Pq) for the cases 
with n = and n = 2. It turns out that after r = 1 the value 
of p is nearly constant forn = 2, but not forn = 0. 

Let us now return to the earlier question of whether a flux 
configuration with zero linking number can have finite spec- 
tral magnetic helicity, i.e. whether H{k) is finite but of oppo- 
site sign at different values of k. The spectra M{k) and H{k) 
are shown in Fig. [9] for the two cases at time r = 5. This fig- 
ure shows that in the configuration with zero linking number 
H{k) is essentially zero for all values of k. This is not the 
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FIG. 9: Comparison of spectra of magnetic energy and magnetic 
helicity in the case with zero Unking number (upper panel) and fi- 
nite finking number (lower panel) at r = 5. Stretcfies with negative 
values of H{k) are shown as dotted lines. 



case and, in hindsight, is hardly expected; see Fig. [9] for the 
spectra of M{k) and k\H{k)\/2iiQ in the two cases at r = 5. 
What might have been expected is a segregation of helicity not 
in the wave-number space, but in the physical space for posi- 
tive and negative values ofy. It is then possible that magnetic 
helicity has been destroyed by locally generated magnetic he- 
licity fluxes between the two domains in > and ^ < 0. 
However, this is not pursued further in this paper. 

In order to understand in more detail the way the energy is 
dissipated, we plot in Fig.[T0lthe evolution of the time deriva- 
tive of the magnetic energy Em = ^ J B'^ dV (upper panel) 

and the kinetic energy Ek = ^ / pU"^ (lower panel). In 
the lower panel we also show the rate of work done by the 
Lorentz force, Wl = J U ■ {J x B) dV, and in the upper 
panel we show the rate of work done against the Lorentz force, 
— VFl- All values are normalized by Euo/Ts, where £^mo is 
the value of Em at r = 0. 

The rates of magnetic and kinetic energy dissipation, eM 
and eK, respectively, can be read off as the difference between 
the two curves in each of the two panels in Fig. [Tol Indeed, 
we have 



Wl - dEu/dt 6M, 



WL^Wc-dEK/dt = eK, 



(12) 



(13) 



where the compressional work term Wc = J p V • U dV is 
found to be negligible in all cases. Looking at Fig. \T0\ we 
can say that at early times (0 < r < 0.7) the magnetic field 
contributes to driving fluid motions (Wl > 0) while at later 
times some of the magnetic energy is replenished by kinetic 
energy (VPl < 0), but since magnetic energy dissipation still 
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FIG. 10: Evolution of the rate of work done against the Lorentz 
force, — Wl, together with d^M / dt (upper panel), as well as the rate 
of work done by the Lorentz force, +Wl, together with d^x/dt 
(lower panel), all normalized in units of Em/Ts, for the case with 
finite linking number. The inset shows — Wl at late times for the 
case with n = (solid line) and n = 2 (dashed line). 



dominates, the magnetic energy is still decaying (dEu/dt < 
0). The maximum dissipation occurs around the time r = 0.7. 
The magnetic energy dissipation is then about twice as large as 
the kinetic energy dissipation. We note that the ratio between 
magnetic and kinetic energy dissipations should also depend 
on the value of the magnetic Prandtl number, PrM = J^/rj, 
which we have chosen here to be unity. In this connection it 
may be interesting to recall that one finds similar ratios of ex 
and eM both for helical and non-helical turbulence |25]. At 
smaller values of PrM the ratio of ex to ex + diminishes 

— 1/2 

like Pr ^ for helical turbulence 12611 . In the present case the 
difference between n = and 2 is, again, small. Only at later 
times there is a small difference in VKl, as is shown in the inset 
of Fig. [Tol It turns out that, for n = 2, VKl is positive while for 
n = its value fluctuates around zero. This suggests that the 
n = 2 configuration is able to sustain fluid motions for longer 
times than the n = configuration. This is perhaps somewhat 
unexpected, because the helical configuration (n = 2) should 
be more nearly force free than the non-helical configuration. 
However, this apparent puzzle is simply explained by the fact 
that the n = 2 configuration has not yet decayed as much as 
the n = configuration has. 
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IV. CONCLUSIONS 

The present work has shown that the rate of magnetic en- 
ergy dissipation is strongly constrained by the presence of 
magnetic heHcity and not by the quahtative degree of knotted- 
ness. In our example of three interlocked flux rings we con- 
sidered two flux chains, where the topology is the same except 
that the relative orientation of the magnetic field is reversed in 
one case. This means that the linking number switches from 2 
to 0, just depending on the sign of the field in one of the rings. 
The resulting decay rates are dramatically different in the two 
cases, and the decay is strongly constrained in the case with 
finite magnetic helicity. 

The present investigations reinforce the importance of con- 
sidering magnetic helicity in studies of reconnection. Recon- 
nection is a subject that was originally considered in two- 
dimensional studies of X-point reconnection [27l |2S|] • Three- 
dimensional reconnection was mainly considered in the last 
20 years. An important aspect is the production of current 
sheets in the course of field line braiding ll29l] . Such current 
sheets are an important contributor to coronal heating [30]. 



The crucial role of magnetic helicity has also been recog- 
nized in several papers However, it remained unclear 
whether the decay of interlocked flux configurations with zero 
helicity might be affected by the degree of tangledness. Our 
present work suggests that a significant amount of dissipation 
should only be expected from tangled magnetic fields that 
have zero or small magnetic helicity, while tangled regions 
with finite magnetic helicity should survive longer and are ex- 
pected to dissipate less efficiently. 
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